Reliable quantum replication at the Heisenberg limit 
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No process in nature can perfectly clone an arbitrary quantum state. But is it possible to engineer 
processes that replicate quantum information with vanishingly small error? Here we demonstrate the 
possibility of probabilistic super-replication phenomena where A'^ equally prepared quantum systems 
are transformed into a much larger number M nearly perfect replicas, with an error that rapidly 
vanishes whenever M is small compared to A*'^ . The quadratic replication rate is the ultimate limit 
imposed by quantum mechanics and is fundamentally linked with the Heisenberg limit of quantum 
metrology. 



No physical process can copy arbitrary quantum states 
on demand [11 [2]: if such a process existed, we could 
build a device that distinguishes quantum states with ar- 
bitrary precision, violating the uncertainty principle and 
enabling faster-than-light communication [3] . Probabilis- 
tic processes like stimulated emission, however, seem to 
evade this restriction. An excited atom interacting with 
a polarized photon is expected to produce sometimes a 
second photon with the same polarization, effectively de- 
livering a perfect clone. True that some other times the 
atom will spontaneously emit a photon of random po- 
larization, but still, when stimulated emission occurs, a 
prefect clone has been produced. If we had on our side 
a quantum version of Maxwell demon, who separates the 
photons produced by stimulated emission from those pro- 
duced by spontaneous emission, we would be able to pro- 
duce any desired number of clones with a non-zero proba- 
bility. Unfortunately, our imaginary helper is not allowed 
by the laws of quantum mechanics: even with small prob- 
ability, it would allow us to violate the no-signalling prin- 
ciple [4|. However, nothing forbids that other probabilis- 
tic processes, akin to spontaneous emission, could boost 
our ability to replicate quantum information beyond any 
previously conceived limit. This possibility raises new 
fundamental questions: Is it possible, for example, to 
engineer a process that duplicates a beam of N equally 
prepared particles, producing a beam of 2N almost per- 
fect clones? What is the ultimate limit to the number of 
replicas that can be produced without significant errors 
by a physical process? 

Here we answer both questions: For a beam of N quan- 
tum particles in an arbitrary state generated by Hamil- 
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tonian evolution, we devise a filtering mechanism that 
allows one not only to duplicate, but also to produce 
an overwhelming number M of nearly perfect clones, 
with a small error that vanishes rapidly whenever M is 
small compared to N'^. We name this new phenomenon 
super-replication and show that it is intrinsically prob- 
abilistic, proving that deterministic processes can only 
produce a negligible number of nearly perfect replicas. 
For example, for 100 linearly polarized photons, super- 
replication allows to produce 1000 replicas with fidelity 
99.9%, whereas the best deterministic process can only 
achieve without filter the fidelity 57%. In addition, we 
show that no physical process, deterministic or not, can 
produce more than O(iV^) replicas: any attempt to repli- 
cate quantum information beyond this limit is doomed to 
produce a joint output that has vanishing fidelity with 
the desired state. To explain the roots of this fundamen- 
tal limitation, we establish a deep link between quantum 
cloning and the precision limits of quantum metrology [S]- 
[TUj, showing that the Heisenberg Limit sets the ultimate 
bound to the replication rate of probabilistic processes, 
while the Standard Quantum Limit sets the correspond- 
ing bound for deterministic processes. 



I. RESULTS 

Quantum metrology limits to information repli- 
cation. The Heisenberg Limit (HL) is the cornerstone 
of quantum metrology [5HTU]. It states that whenever 
one tries to estimate the value of a parameter charac- 
terizing a physical process, the variance of the estimate 
will decrease at most as 1/iV^, where N is the number 
of times that the process is probed. The scaling 1/iV^ 
can be achieved, provided that the probes are prepared 
in a suitable entangled state. However, if the probes are 
not entangled, then the variance decreases only as 1/7V — 



known as the Standard Quantum Limit (SQL), which is 
the scaUng of a classical statistics over repeated experi- 
ments. 

At first sight, the Heiscnbcrg Limit seems to have noth- 
ing to do with quantum cloning: the former limits the 
quality of estimation using entangled states, the latter 
aims at producing copies from a product state in input. 
However, we will now show that the HL sets the ultimate 
upper bound on the performance of cloning. We consider 
the setting where N identical systems are given, and each 
system is prepared in the same state \^pg) = e^^^\ip), for 
some Hamiltonian H and for some unknown parameter 
6' € M. For example, 6 could be the polarization angle 
of a linearly polarized photon, or the time parameter in 
a dynamics governed by the Schrodinger's equation. To 
link quantum cloning and quantum metrology, we adopt 
an asymptotic approach inspired by information theory: 
we consider cloning processes that transform N input 
copies into M = cN°' approximate copies, where c > 
and a > 1 are suitable constants. We call the constant a 
the replication rate and say that a process achieves reli- 
able replication if the error vanishes in the limit of large 
N. 

Now the question is: how many copies of the state lipg) 
can be produced by a reliable replication process? The 
answer, argued for in the Methods section, is given by 
the following: 



1. Heisenberg Limit for Quantum Replication: 

No physical process can reliably replicate the quan- 
tum states {e^^^\ip),6 G M.} at a rate larger than 
2. 

The Heisenberg Limit leaves lots of room for replicat- 
ing quantum information: for every rate 1 < a < 2 one 
has the chance not only to duplicate the input copies, 
but also to produce an overwhelming number of replicas. 
In the next paragraph, we will exhibit explicit schemes 
that achieve replication at every rate allowed by the HL. 
These schemes are necessarily probabilistic, as a conse- 
quence of the Standard Quantum Limit: in the Methods 
section we show that the SQL implies that every deter- 
ministic process is bound to produce a number of replicas 
that is negligible with respect to N: 

2 Standard Quantum Limit for Quantum 
Replication: No deterministic process can reliably 
replicate the states {e^^^\tp), 9 €R} at a rate larger 
than 1. 

In other words, without a probabilistic filter, one can only 
embezzle from nature a negligible number of extra-copies. 
The difference between probabilistic and deterministic 
processes does not consist only in the number of replicas 
that can be produced reliably, but also in the speed at 
which the replication error vanishes: as proven in the 
Methods section, the error cannot vanish faster than (the 
inverse of) a polynomial of degree smaller than 4 in the 



deterministic case, whereas in the probabilistic case the 
error can vanish faster than every polynomial. 

Probabilistic super-replication. Here we show that 
the states \ipe) — e*''^|^) can be reliably replicated at 
any rate allowed by the HL. When the rate is larger than 
(or equal to) 1, we call the process super-replication to 
emphasize the fact that it beats the SQL. In addition, 
we pin down the probabilistic nature of super-replication, 
showing a tight relation between the replication rate and 
the probability that replication takes place successfully. 
The key idea to achieve super-replication is to devise a 
quantum filter that modulates the Fourier amplitudes of 
the wavcfunction, in a way that enhances the cloning per- 
formances to the maximum rate allowed by quantum me- 
chanics. To move to the Fourier picture, we express the 
input states as {tpe) = Z]_EGSpec(ff) v^^'^"^!-^)' where 
Spec{H) denotes the spectrum of the Hamiltonian, \E) 
is an eigenvector of H with eigenvalue E, and pe is the 
probability that a measurement of the energy gives out- 
come E. Without loss of generality, we assume that the 
average energy is zero, a condition that can always be 
enforced by shifting all the eigenvalues by a suitable con- 
stant. When N identical copies are given, the joint state 
can be expressed as 
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e''^y/f^\N,E), (1) 



where H^^'' — X]i=i ^i i^ the total Hamiltonian (Hi de- 
noting the one-body Hamiltonian H for the z-th system) , 
I A^, E) is an eigenvector of H^^> for the eigenvalue E, and 
Pn,e is the probability that a measurement of the total 
energy gives outcome E. 

The ideal goal of cloning is to produce M perfect copies 
in the state \ipg)^'^ ■ The filter is described by two mea- 
surement operators My^s and Mno, acting on the space 
of N copies. When the N input systems emerge from 
the filter they are in the state Myes\ipe)^^ , up to a nor- 
malization factor. Now, by suitably choosing the oper- 
ator Myes, we can modulate the Fourier amplitudes of 
the input state in any way we like. One natural choice 
would be to replace N with M in the probability dis- 
tribution pn,e- Since, however, the spectrum of H^-^' 
may be not be contained in the spectrum of H^-'^\ we 
may need also to shift the values of the energy by a small 
amount Eq, which tends to zero for large M (see the 
Supplemental Material for details) . With this choice, the 
state after the filter is projected to the entangled state 

l*e^) « EBgSpec(ff««) e^'^'y^PM^E+EoW, E) . 

The state |^^) will now act as a quantum program, 
containing the instructions that will be used by a de- 
terministic quantum device to generate M approximate 
clones. For this purpose, we consider a quantum device 
that coherently transforms each basis state \N,E) into 
the corresponding state |M, E + Eq), thus producing the 
state |$f ) (X EiJeSpec(//(")) e'''' ^Pm,e+e,\M, E + E^). 
It is not hard to see that the state |<I>g^) has a very high 



overlap with the desired state \ipe)'^^ for every rephca- 
tion rate aUowed by the Heisenberg hmit. For large iV, 
using the central limit theorem one can show that the 
fidelity satisfies the following lower bound 
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where r/ is a suitable constant (see Supplemental Ma- 
terial). The bound shows that the fidelity rapidly ap- 
proaches 1 whenever M is of order TV" with a smaller 
than 2, that is, whenever the replication rate satisfies 
the Heisenberg Limit. A strong converse can be proven 
too: every physical process replicating quantum states at 
a rate higher than 2 must have vanishing fidelity in the 
limit of large N . 

Thanks to the use of our filter, we have been able to 
embezzle from nature a large number of replicas. The 
improvement is striking if one compares it with the per- 
formances of standard cloning processes. Let us illustrate 
this fact for the replication of linear polarization states 
l^e) = cos9\V) + sa\9\H). In this case, the best de- 
terministic process is the phase-covariant doner of Ref. 

[TT] . and its fidelity is F'^^^,^ « "^m+n '^"^ ^^^ asymp- 
totic limit of large N and M . In agreement with our 
Standard Quantum Limit, the fidelity vanishes when- 
ever M is of order N^'^'^ , e > 0: As we proved in the 
general case, deterministic processes can only produce 
a number of replicas that is a negligible fraction of the 
number of input photons. For example, starting 100 
from input photons, a deterministic process can produce 
only a few extra-copies, while using our filter one can 
produce up to thousands nearly perfect replicas. For 
N = 100 and M — 1000, cloning with filter achieves 
fidelity Fprob — 0.9986, while without filter the fidelity 
is Fdet = 0.5739. The different features of replication 
processes at the Heisenberg limit and at the standard 
quantum limit are illustrated in Fig. [T] A plot compar- 
ing the performances of replication with and without the 
filter in the case iV = 20 is presented in Fig. [2] What 
makes the improvement even more dramatic is that the 
quality of the replicas is measured by the global fidelity 
between the output state and the desired joint state of 
M perfect copies, which is much smaller than the overlap 
that each single copy has with the state |'0e)- 

Maximizing the probability of successful repli- 
cation. The dramatic enhancement leading to super- 
replication comes at a price: the probability that the 
input systems pass the filter has to decay with N. For 
our filter, the probability decreases exponentially fast as 
Pyes[N — > M] < e~^^ for a suitable constant k > Q. 
One can ask whether there are ways to achieve super- 
replication with a larger success probability. It turns out 
that the answer is yes, and that larger success proba- 
bilities can be obtained through a minimal modification 
of our filter. In particular, one can generate M = cN 
replicas with a probability Pyes [N — )■ M] going to zero as 
1/N for every desired 5 > 0. The price to pay, however, 
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FIG. 1. More and better: the advantages of proba- 
bilistic super-rephcation (a) Every deterministic process 
is constrained by the Standard Quantum Limit: the number 
of clones that can be produced reliably from A*' input copies 
is of order A^ and the fidelity of the clones cannot approach 
1 faster than the inverse of a polynomial of degree 4. (b) 
Super-replication at the Heisenberg limit: the cloning perfor- 
mances can be dramatically increased by a probabilistic filter, 
depicted here as a "quantum Maxwell demon" that separates 
two branches of the wavefunction. In the successful branch, 
any number of clones of order A'^'^"'^ can be produced with 
fidelity approaching 1 faster than any polynomial. The prob- 
abilistic filter allows one to generate more clones and with a 
better quality. 




FIG. 2. Rephcation of linearly polarized photons: 
with and without filter. The fidelity of the best cloning 
processes for linearly polarized photons \tpg) = cosO\V) + 
sm9\H) is plotted here for A'' = 20 input copies, with the 
number M of output clones ranging from 20 to 400. The red 
(green) line refers to optimal cloning with (without) filter. 



is that the fidehty will go to 1 only as 1/{N^ hvN). More 
generally, in the Supplemental Material we show that 
every replication process with non-vanishing probability 
must necessarily have vanishing fidelity. This provides a 
strong converse of the Standard Quantum Limit: every 
deterministic process producing replicas at a rate higher 
than the SQL must have vanishing fidelity. 

For super-replication processes with rate a — l-l-e, e > 
0, we show how the success probability can be increased 
to Pyes[N -^ M] > e^^ for every desired exponent 5 > 
e. However, no further improvement is possible below the 
critical value 6c = e: any process with success probability 
scaling as Pyes[N -^ M] — ae~^^ , 5 < e and for a,b > 
must have vanishing fidelity in the asymptotic limit. 

Delayed choice between deterministic and prob- 
abilistic strategies. If information is money, an easy 
way to picture the difference between deterministic and 
probabilistic replication is as two different ways to make 
investments. Deterministic processes are like depositing 
money in a bank with low interest (one gets a negligible 
amount of extra replicas), whereas probabilistic processes 
are like buying tickets of a lottery with large jackpot (the 
prize is big, but the chance to win is small). Depending 
on the situation, one or the other approach may be more 
convenient. Interestingly, a bonus of our results is that 
one can always delay the choice whether or not to use 
the filter. Indeed, one can first implement the determin- 
istic transformation \N,E) i— )■ \AI, E + Eq) on the energy 
eigenstates and postpone the modulation of the proba- 
bility distribution to the very end. At this stage, if the 
modulation is omitted, one has a nearly optimal deter- 
ministic process for replication at the SQL. Instead, if 
the modulation is performed, one can probabilistically 
enhance the replication performance to the Heisenberg- 
limit scaling. 

The many-worlds fairness of quantum mechan- 
ics. Super-replication can be achieved not only for one- 



parameter families of states of the form \ipg 
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but also for different manifolds of states, including the 
manifold of all maximally entangled states of two identi- 
cal systems. However, super-replication is not a generic 
feature. For example, if one tries to copy an arbitrary — 
as opposed to linear — polarization state, then no filter 
is going to help: the performances with filter are equal 
to the performances without filter. What is the origin 
of such a difference between linear and arbitrary polar- 
ization? As we show in the Supplemental Material, the 
explanation comes from the geometry of quantum states: 
states on a circle can be replicated at the Heisenberg 
limit, but states on the sphere cannot. We also show 
that such a limitation holds for quantum systems of ar- 
bitrary dimension: no probabilistic filter can make an 
arbitrary quantum state more copiable. This property 
is quite compelling when considered from the angle of 
the many-worlds interpretation of quantum mechanics, 
because it states that no-branch of the wavefunction of 
the universe offers an advantage over the others when it 



comes to replicating the information contained in a com- 
pletely unknown state. Regarding different branches as 
"different worlds" , we can formulate this as a fundamen- 
tal principle, which we name many-worlds fairness, stat- 
ing that the capability to replicate arbitrary information 
is the same in all possible worlds. Many world fairness 
rules out quantum mechanics on real Hilbert spaces |12j , 
an alternative physical theory where photons can have 
only linear polarizations. Thanks to this observation, one 
can provide a new answer to the old question "Why are 
physical systems described by complex, — instead of real 
Hilbert spaces?" . Traditionally, the standard answer has 
been to invoke local tomography, namely the property 
that one can completely identify a mixed state from the 
statistics of local measurement on the components. How- 
ever, one may consider it as an ad hoc requirement .15] , 
and, in fact, there are even reasons to prefer real quan- 
tum mechanics to its complex version, as it was recently 
pointed out by Wotters JJj. Balancing this fact, many 
world fairness offers a new reason (other that the usual 
local tomography) to prefer complex quantum mechanics 
to its real relative. 



II. DISCUSSION 

Towards experimental implementation. Super- 

replication at the Heisenberg limit has been introduced 
here from a theoretical point of view. But is it possible 
to implement it experimentally? Luckily, very recently 
there have been important experimental breakthroughs 
on the closely related topic of probabilistic amplifica- 
tion of coherent states of light [TSKIQ] . Although the 
translation to our case is not immediate, we suggest that 
super-replication of linearly polarized photons could be 
achieved through a suitable sequence of amplitude damp- 
ing channels, which for the polarization play the role of 
the photon subtraction for coherent states. An alterna- 
tive approach would be to first encode the states of the in- 
put copies into coherent states via matter-light teleporta- 
tion [50] , amplify the latter, and then teleport back. This 
would provide a new application to the existing experi- 
mental schemes of coherent state amplification, making 
them the building block for the replication of quantum 
information at the Heisenberg limit. 

Relation with previous works on probabilistic cloning. 
The idea that in some situations even non-orthogonal 
states can be copied perfectly by a probabilistic device 
was first introduced by Duan and Guo TlJ, who showed 
that a set of pure quantum states can be copied per- 
fectly if and only if they are linearly independent. This 
means that for a single photon one can perfectly copy two 
polarization states, while for N photons one can copy at 
most A^+l states. Subsequent works [221I13] showed that 
probabilistic devices can improve the performances of ap- 
proximate cloning, leading to nearly ideal performances 
in the case of coherent states with equal amplitude. The 
nearly perfect cloning of coherent states may lead one 



to believe that there is no hmit to the amount of clones 
that can be produced probabilistically. Contrarily to this 
intuition, we have proven here that the Heisenberg limit 
sets the ultimate bound M = 0{N'^) to the number of 
clones that can be produced reliably. Moreover, we have 
shown that for arbitrary one-parameter families of states, 
any replication rate allowed by the Heisenberg Limit can 
be achieved, and that the process that achieves it must 
be probabilistic, because deterministic processes are con- 
strained by the Standard Quantum Limit M = 0{N). 



Relation with probabilistic quantum metrology. Recently 
it has been observed that the precision in the estimation 
of a parameter encoded in a set of quantum states can 
be improved using probabilistic schemes |24( 125) where, 
depending on the outcome of the filter, one can decide 
to abstain from estimating the parameter. The mecha- 
nism of abstention can sometimes boost the precision of 
estimation from the standard quantum limit 1/iV to the 
Heisenberg limit 1/iV^. The connection between these 
results and the super-replication phenomenon is made 
clear by our approach to the Heisenberg limit: the fact 
that a probabilistic filter can improve estimation until 
the Heisenberg limit implies that N'^, rather than N, is 
the upper bound for the replication rate of the states 
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\ip). However, there is an important catch: 



the filter that achieves replication at the Heisenberg limit 
is not the same filter that achieves estimation at the 
Heisenberg limit. In fact, if we were to use the same 
filter needed in Refs. [241 ES], we would not be able to 
reduce the error down to zero. 



Conclusions. We discovered a previously unexplored con- 
nection between quantum cloning and quantum metrol- 
ogy, showing the Standard Quantum Limit and the 
Heisenberg Limit for estimation can be translated into 
limits for the reliable replication of quantum information. 
The possibility of producing a large number of almost 
perfect replicas, opened by the Heisenberg Limit, has 
been demonstrated for arbitrary one-parameter families 
of states generated by a Hamiltonian H, showing that, 
with the help of a probabilistic filter, one can embezzle 
from nature a large number number of extra clones, as 
long as the number is small compared to A'^^, the square 
of the number of input copies. We then showed that the 
phenomenon of super-replication is not generic: In par- 
ticular, there is no way to improve the performances of 
cloning of an arbitrary polarization state. The fact that 
the geometry of quantum states forbids any advantage in 
the replication of an arbitrary state based on the prob- 
abilistic branching of the wavefunction is a fundamental 
property of quantum mechanics, that we call many-world 
fairness. It is our hope that our work will stimulate the 
debate on the foundations of quantum mechanics and 
trigger new advancements on the experimental realiza- 
tion of probabilistic quantum information-processing de- 
vices. 



III. METHODS 

Decomposition of quantum instruments. To de- 
scribe the most general probabilistic processes allowed by 
quantum mechanics, we use the framework of quantum 
instruments. A quantum instrument with input (output) 
Hilbert space J^m {J^out) is a collection of completely 
positive, trace non-increasing linear maps {"PijigXi each 
map transforming density matrices on J^„ into sub- 
normalized density matrices on J^out- If the input system 
is prepared in the density matrix p, the probability that 
the i-th process takes place is p{i\p) = Ti['Pi{p)]. Condi- 
tionally to the occurrence of outcome i, the state of the 
output system is p'^ = Vi{p)/ Tr[Pi{p)]. We note an ele- 
mentary fact about quantum instruments: any quantum, 
instrument can be decomposed into a pure measurement 
on the input .system followed by quantum channel depend- 
ing on the outcome. By "pure measurement on the input 
system" , we mean a quantum instrument {A^ijigx of the 
special form Aii{p) ~ MipAlJ for some measurement op- 
erator Mi on Jifin- By "quantum channel", we mean a 
completely positive trace-preserving map Ci, transform- 
ing density matrices on J^„ into (normalized) density 
matrices on J^out- Hence, our claim is that every instru- 
ment {T^ijigx can be decomposed as 



V^ = C,M^ 



VieX. 



(3) 



The proof is provided in the Supplemental Material. The 
decomposition of Eq. ([3]) is important because it states 
that any probabilistic process in quantum theory can be 
decomposed into a measurement, which modulates the 
amplitudes of the input state in a suitable basis, followed 
by a deterministic transformation depending on the out- 
come. In our case we consider measurements with two 
outcomes, referred to as filters. The filter induces a bi- 
furcation of the wave-function and selects one particular 
branch, corresponding to one outcome with measurement 
operator Aly^s- 

Quantum metrology bounds. Here we sketch the 
proof that deterministic processes cannot reliably pro- 
duce more than 0(7V)replicas (SQL), whereas probabilis- 
tic processes cannot produce more than 0{N'^) replicas 
(HL). We present the main ideas, while the full detail 
can be found in the supplementary material. Consider 
an exponential family of quantum states, of the form 
(IV'e) — e^^^\ip)}. When N identical copies are avail- 
able, the quantum Cramer-Rao bound [7] states that 
the variance in the estimation of 9 is lower bounded 
by Vm{0) > ^, for a constant c which can be set to 
c — I with a suitable change of parametrization. In 
this case, it is known that the bound can be achieved 
for large N. Suppose now that a deterministic pro- 
cess C produces M = N + N", with an error vanish- 
ing in the limit of large N. To measure the error, we 
use the trace distance between the output of the process 
and the desired state of M perfect copies. For simplic- 
ity, we assume here that the error can be bounded as 



\\C{\i;e){iPe\'^^) - miM'^^h < ^-^ for some /3 > 1 
(this requirement will be lifted in the Supplemental Ma- 
terial, where we prove strong converses to the HL and 
SQL). We now use fact that the trace distance gives a 
bound on the total variation distance between the proba- 
bility distributions arising from quantum measurements. 
This implies that the variance in the estimation of 9 
from the clones, denoted by Vout{6)^ can be bounded 
as Vout{0) < Vm{0) + 7M~^, for a suitable constant 
7 > 0. However, applying the quantum channel C can- 
not reduce the variance of the optimal estimation strat- 
egy, denoted by V^'"(6'). Therefore, we have the bound 
¥^'"(6*) < Vout{d) < VM{9)+-fM-^. Choosing the best 
estimation strategy, we then obtain 



VS''''{9)<V^f''{0)+jM- 



(4) 



For large N, we can replace the variances with the right 
hand side of the Cramer-Rao bound, thus obtaining 



A/N < 4/(7V + N") + -f{N + N")-^. 



(5) 



Clearly, this implies a < 1, because otherwise we would 
have a contradiction for large N. This establishes the 
standard quantum limit for reliable cloning. Moreover, 
we can also show that the reliability exponent /3 cannot 
be larger than 2. Indeed, by Taylor-expanding the r.h.s. 
of Eq. (Is]) we obtain the condition [3 < 2 — a, which 
implies /? < 2. This means that no deterministic cloning 



process can produce copies with a trace distance vanish- 
ing faster than iV~^. Equivalently, the fidelity cannot 
approach 1 faster than N^"^. The difference with prob- 
abilistic processes is striking: remember that our filter 
allows to clone polarization states with an error that van- 
ishes faster than any polynomial. 

The derivation of the Heisenberg limit for reliable 
cloning follows the same lines sketched above. In the 
probabilistic case, the input state of the quantum channel 
Cyes is not the product state \tpg)'^^ , but instead the en- 



tangled state \^g) = MyeslV'e 



\tg>N , 



\My,,\^e)^^\\. Tak- 



ing the average over all possible values of 9, the Heisen- 
berg limit implies that the variance will be bounded 
as Vin.ave > c/N^, for a suitable constant c > (see 
the Supplemental Material for the proof). Choosing a 
cloning rate a > 1 and a reliability exponent /3 > 1 
and substituting the values in the quantum Cramer-Rao 
bound, Eq. H then becomes 



cN-^ <4:N-" + 0{N-°'), 



(6) 



which implies a < 2. Note that here there is no bound on 
the reliability rate (3. This feature enables the exponen- 
tial scaling of the fidelity demonstrated by our strategies. 
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Appendix A: Quantum Cramer-Rao bound for probabilistic metrology 

Here we establish the Heisenberg limit for the probabihstic estimation of the parameter 6 characterizing the state 
li'e) = e***^!^); showing that the variance can decrease at most as iV~^ with the number iV of input copies. In order 
to prove this result, we have first to estabUsh a suitable Cramer-Rao bound for probabilistic quantum metrology. We 
will first by recalling some known facts about the Cramer-Rao bound in the usual deterministic scenario, and then 
use these facts to treat the probabilistic case. 

In quantum metrology, the variance in estimation of parameter 6 from a pure state j^e) is lower bounded by the 
Cramer-Rao bound 



v{e)^ 



1 



FQil^e)) 



Fgil^e)) ■■= ^AHf 



(Al) 



where Fqd'^g)) called quantum Fisher information, H is the Hamiltonian that governs the evolution of l^e), namely 
i^l^e) = -ff|*e) and {AH)^ := (^ejiJ^j^g) _ (*e|i/|*e)2 is the variance oi H. 

For A'' copies of a pure state \ipe) — e'^^j-fA), the Hamiltonian is given by iJ^^^ — 'J2i=i^jj where Hj denotes 
the Hamiltonian H acting on the Hilbert space of the j-th copy. Hence, the quantum Fisher information is given by 
Fgdipg)'^^) = 4N{AH)^. Inserting this value in the Cramer-Rao bound one obtains the standard quantum limit, 
stating that the minimal variance for estimation using the product probe state l^e)®^ scales as 0{N~^). 

We now consider the probabilistic case, where we want to estimate the parameter 9 from the state j^g) using a 
filter. The joint state of the systems emerging from the filter is given by 



l$fl 
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To compute the Fisher information, we need to compute the variance of the Hamiltonian Kg defined by i^l'&e 



Kg\^g). Note that, in principle, Kg 
calculation we obtain 



can depend on 9, even if the original Hamiltonian H does not. By explicit 



Kgl^g) = 



My^MM-^'^ 
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The quantum Fisher information at 9 is then given by 



FQ{\<^g)) = {>^e\KgKg\>^^) ~ |($e|i^e|$ 



(A4) 



Now, in order to minimize the variance one has to maximize the Fisher information over all possible filter operators 
Myeg. In this problem, there is an important caveat: as a matter of fact, for every given value of 9 there one can choose 
a suitable My^g that makes the quantum Fisher information as large as desired. This can be easily seen by choosing 
a parametrization where {'^g\H\'^g) = and defining My^s.e = e|^e)(*e| — l^e")(*e"l' where |\l/g-) is the normalized 
vector proportional to H\^g). This choice gives a Fisher information scali ng as 1/e^, which can be made arbitrarily 
large by choosing e to be small. Hence, the usual Cramer- Rao bound (Al) would imply that probabilistic metrology 
can make the variance as small as desired. However, this is an artifact of the fact that the usual Cramer-Rao bound 



is a local bound: a probabilistic estimation scheme can have an unlimited precision only at a particular value of 9, 
which is the unknown parameter. Since 9 is unknown, instead of minimizing the variance at a particular value of 9, 
it is more sensible to minimize the average variance. 



Vave ■■= I TT P(^lyes) V (9) 



where Q is the set of all possible values of 9, p(9\yes) is the conditional probability distribution 
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and p{9) is the prior distribution of 9. For simplicity, in the following we will consider the case of periodic evolutions, 
with O — [0, 27r), in which case we choose the prior to be uniform (i.e. p{9) — 1 for all 9). In the end of the section, 
we will discuss how to modify the argument in the non-periodic case = JR. 

It is easy to see that the average variance is lower bound ed b y the inverse of the averge quantum Fisher information 



Fave ■= jQd9/{27r) p{9\yes) Fgd^g)). Indeed, using Eq. (Al) we have the 
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which implies the Cramer-Rao type bound Vave > jr-^ — • Minimizing the variance over all filter operators Myes we 



then obtain 



fQ.a^ 
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(A5) 



This form of the Cramer-Rao bound is the most suitable to study probabilistic metrology. We will now use it to 



IV'6 
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\ipd) = e \ip), the average quantum Fisher 



show that for a family of states of the form (A2), with j^I^e 

information scales as 0{N'^). For simplicity, we will consider the case of phase estimation, where H has integer 
spectrum and 9 is an angle 9 e [0, 27r), which is the case considered in the main text. Discarding the second term in 
the r.h.s. of Eq. (A4) and using Eq. (A3) we obtain the upper bound 



FQ{\^e))^^{^e\KlKg\^e) 
= ^Ke\^e)\? 



<4 



^4 



M,e.7?(^)M-l|<i>,) 



+ 



(*e|[i?W,Mt,,Af,e.]|* 



yiyes\'i!e) 



\My,sH^^)\^e)\\ , \{-^e\Ml,My,MW\ii,,) 



^4 



\\Myes\^e)\\ 

'2||M,,,i?W|*,)|| 



\\Myes\^eW 



\\Myes\^e)\\ 
(*e| MlesMy^,\^g) 



Taking the average over 9 we then have 
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Noting that Myes commutes with iJ'-^' , we then obtain 

FnZe < max ^ ye^^^ ^ 

16 (*! MlsH^^'^H^^^ Myes I*) 
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where H-ff'-^'Hco denotes the operator norm of H^^', equal to the maximum of the modulus of its eigenvalues. Since 
jj{N) _ ^ - -^ Hj we have that ||i?^^-'||oo scales linearly with N, thus giving the bound Fq'^^ < 0{N'^), and therefore 

K'S^" > 0{N-^). 

In the non-periodic case 9 E R, all the arguments can be repeated by choosing a Gaussian prior px{9) = 
■\/27r/cr2e~^ /'-^'^ \ the only difference being that now the operator Myes does not commute with H^-'^K However, 
we can always replace Myes by its decohered version, where all the off-diagonal elements in the eigenbasis of H^^^ 
are replaced by zero. Since these off-diagonal elements are exponentially small in the variance cr^, the corresponding 
correction can be made arbitrarily small, leading to the bound Fq^^ < 16\\H'-^^\\-^ + e with arbitrary small e. 

Appendix B: Decomposition of quantum instruments 

Here we prove that any quantum instrument {T'iliex can be decomposed into a pure measurement on the input 
system, followed by a deterministic process conditional to the outcome of the measurement. 

The proof is simple: let Vj be the adjoint of the map Vi, defined by Tt[PJ{A)p] = Tr[AVi{p)] for every operator 

A and for every density matrix p. Define the operator Mi :— \/'Pj{F} and the map Ci by the relation Ci{p) :~ 



Vi{M^^pMj^^), where M^^ is the inverse of Mi on its support. If Mi is invcrtiblc, then Ci := Ci will be trace- 
preserving (and hence a quantum channel). In this case, by definition we will have Vi — CiAii. If Mi is not invertible, 
we can define the trace-preserving map Ci{p) := Ci{p) + QopQoj where Qo is the projector on the kernel of Mi. Again, 
the definition implies Vi — CiAii- Hence, we established the decomposition Vi = CiAii, Vi G X. 

Appendix C: General facts about quantum cloning 

Here we give the general formulation of the problem of optimal cloning in the presence of symmetry, taking into 
account arbitrary probabilistic processes. The general facts introduced here will be specialized to the case of cloning 
of states of the form \^pg) = e'^^^\ip) in the next section. 

1. Probabilistic cloning in the presence of symmetries 

We consider probabilistic cloning processes designed to produce approximate copies of a pure state of the form 
IV'g) := Uglip), where \ip) e J^ is a fixed unit vector and {Ug} is a unitary representation of a group G on the 
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Hilbert space of the system. A probabilistic cloning process is described by a quantum instrument {Pyes,'Pno} 
consisting of two quantum operations (completely positive trace non-increasing maps) corresponding to the successful 
and unsuccessful instances of the process, respectively. For a cloning machine designed to transform N input copies 
into M approximate copies, the probability of success is given by p{yes\g) := T'i:{'Pyes[{\'>Pg) {''I'gl)'^^]} and the state 
of the M output systems is p'g := 7'ye.[(IV'<,>(V'ffl)®'^]/Tr{Pye«[(|^g)(V'sl)®'^]}- 

When the group G is compact, it is natural to assume that the prior probability of the states Itpg) is uniform and 
given by the haar measure g. In this case, the average fidelity of the output of the device with the desired state 
|^^^(giM^ conditionally to the successful realization of the process, is given by 






F ^ gp{g\yes) {^Pgr' p'Mr'' = ' ^7' ^ ^ .?/ w . Jn ■ (CI) 



From the form of Eq. ( CI ) is is clear that the search for the optimal quantum operation Vyes can be restricted without 
loss of generahty to the set of maps satisfying the covariance condition VyesUf^ = l^f^'Pyes.^g S G, where Ug is 
the unitary channel defined by Ug{p) — UgpU^. When the group G is non-compact, restricting to covariant processes 
means requiring that both the probability of success and the fidelity are independent of the unknown group element 
9- 

Decomposing the covariant quantum operation Vyes a-s Vyes — CyesM.yes where M.yes{p) — ^^yesP^'^les i^ a pure 
measurement and Cyes is a quantum channel (completely positive trace-preserving map) we have that the latter must 
satisfy the symmetry conditions 

[Myes,Uf'']=Q ygeG (C2) 

CyesUf'' = Uf^Cyes V.g e G. (C3) 



Using the decomposition Vyes = CyesM.yes and the above symmetry conditions, one can transform Eq. (CI I into 
the form 

F = {ijf^'Cyes{\^''){^''\m^^'. (C4) 

where \^^) is the conditional state 

^ ._ Myes\i^r'^ .C5) 

I* ^-WMyesm^^W ^^'^ 



Note also that the symmetry condition of Eq. ( C2 ) implies that the probability that the input system pass the 
filter is independent of g and is given by pyes = ||Myes|V')^^lP- 

When {t/g, g € G} are all the unitary operators on the Hilbert space of the system, any operator Myes that 
commutes with all unitaries will act on \ipg)^'^ as a multiple of the identity, due to the Schur's lemma. Hence, 
the action of the fiher will not change the input state: |$^) = Myes|V'ff)®^/IIIV'ff)®^ll = IV'g)'^^ for every value 
of g. We conclude that in this case the symmetry of the set of states inhibites super-replication, by enforcing the 
standard quantum limit. The same conclusion is reached if {C/g,g € G} is an irreducible representation and if |^) is 
a lowest weight- vector, because also in this case the operator Myes acts as a multiple of the identity on the states 
jf/'g)®^, inhibiting super-replication. For example, super-replication is impossible for the arbitrary coherent states 

|a) — e^l"l'/^5I]i^o'^"/"^^l"') *^^ *^® electromagnetic field (as opposed to the case of coherent states with fixed 
amplitude) . 

Appendix D: Super-replication of Hamiltonian-generated quantum states 

Consider the case the group is either C/(l) ~ [0, 27r), in the case of periodic evolutions, or M, in the case of aperiodic 
evoultions. The representation of interest is Ue — e'^^, and the input state is \tp) = X^BeSpecfi?) %/Ps|^), where 
Spec(iJ) denotes the spectrum of H and \E) S Jf is a normalized eigenvector corresponding to the eigenvalue E. We 
will assume without loss of generality that the vectors {\E) \ E E Spec(iJ)} span J^f and that {ip\H\il}) — 0. When N 
copies are given, we write IVj)**^ — X^BeSpecCZ/f")) \/Pn,e\N, E), where \N, E) £ J^®^ is a normalized eigenvector of 
H^^i corresponding to the eigenvalue E. By the central limit theorem, for large N the probability distribution pn.e 
converges in expectation values to the Gaussian gN.E = e~^ /{2N{h )) j ^ 21: N [W^) . 
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1. strong converse of the quantum metrology limits 

Here we consider super- replication processes producing M = cN^^'^ , e > 0, c > output copies. We determine the 
ultimate limits to super-replication, showing that 

1. Every super-replication process with success probability py^s > ae~^^ for some a,b > and x < e must have 
vanishing fidelity. In particular, this implies that every deterministic process producing more than 0{N) output 
copies must have vanishing fidelity. 

2. Every probabilistic process producing more than cN'^ output copies for some c > must have vanishing fidelity. 

The above results provide a strong converse of the standard quantum limit and of the Heisenberg limit, respectively. 
Let us introduce the Cho i op erator of the channel Cyes, denoted by Cyes- In terms of the Choi operator, the 
symmetry condition of Eq. (C3) becomes [Cyes,Uf^ (^ t/®g^] = 0, V0 £ [0, 27r). Hence, we have that Cyes can be 
written in the block diagonal torm Cyes = ® pc eciHC^'i-HC^)] ^a" with each block in the form of 



Cf.= J2 i(^f^]E,E' \M,E' + ^i){M,E + ^i\®\N,E){N,E'\, (Dl) 

where the summation ^ is restricted to the values of E (£") such that E + ii (£" + /i) belongs to Spec{H^^^^). 
The trace-preserving condition on Cyes can be re-written as 



J2 [C^]e,e = 1 Vi? e Spec (i/W) , 



(D2) 



where the summation is restricted to the values of /z such that E + /j. belongs to Spec(iJ). The fidelity in Eq. (C4) 
becomes 

F={^/j\'^^''{^^\Cyes\^r^\^''}, 



to be maximized over the set of positive operators Cyes satisfying the constraints of Eqs. (Dl ) and (D2 1. Expanding 
the states Ith)'^^^ and 1*^) as 



EeSpec^Hl"')) 

l*^>= E '^e\N,E), 

£;eSpec(H(")) 

(for some suitable coefficients {vP^} satisfying X]_EeSpecf _?/(«)) I^bI^ — 1)' ^^ then obtain 



(D3) 
(D4) 



^ = E E [Cf.]E,E'VpiF:^^%*vpif:¥T;i^§> 

P \e,E' 



(D5) 



The positivity of the operator Cyes (and hence of the matrices C^) implies the inequality |[Cp]£;_£;'| ^ \/<^e'^e'^ 
Cfi •= [Cfi]E,E, which leads to the bound 
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fJ. L E 



-\ 2 
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Picking < 5 < e and defining Es = 2^/W+^{lP) we get 
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For large N, the first term in Eq. (D7) is vanishes whenever S < e: indeed, one has 



2 



M V|B|<B. V ^2^- / ^^°° M V|£|<£. / \\E\<Es Py^' 

^„'™EfE»«]fE ^"'^ 

/^ \\E\<Es I \\E\<E, ^V"" 

\,E\<Es 

where the last equality is due to the fact that 'Pm,e converges to a Gaussian with variance of order viV^, and the sum 
is over a set of negligible size ViV* . 

If e > 1, this is enough to prove that the fidelity vanishes, no matter how small the probability is: indeed, choosing 
1 < (5 < e, we have that \E\ < Es for every E £ Spec(ij(^^) and therefore, we do not have the second term in Eq. 



(D7). Hence, we have proved a strong converse of the Heisenberg limit: every process producing more than 0{N'^) 



output copies must have vanishing fidelity. 



If e < 1, then the second term in Eq. (D7) can be upper bounded as 

2 



gN,Es sr^ I \r^ ic%PM,E+ti 



iim y^ y^ TTE\/c'^PM.E+tiPN.E < iim — y^ y^ ,, 

M \\E\>Es I ^y"" M \\E\>Es V ^y^ 



JV->-00 Py^g 



M \E>Es / \E>Es 



9N,Es (N+\S^ec{H)\~l 



- N^oo py^, \ |Spec(iJ)|-l 

1 ,. e 

< — , lim 



^27r(iJ2) W^oo py^ 



„-N 



Now, \i Pyes > ae ^^ for some a, 6 > and x < e, then choosing 5 such that x < 5 < ewe have lim7v_>.oo ^ = 0, 

implying that the fidelity vanishes. This proves a strong converse of the standard quantum limit: every deterministic 
process producing more than 0{N) output copies must have vanishing fidelity. 

2. Performances of the filter 

Here we evaluate the fidelity and the probability of success of the super-replication strategy proposed in the main 
text. Let us first recall the ingredients: The probabilistic filter is given by 

My,,= Yl ^e\N,E){N,E\ tte = 7,/^^^^^±^, (D8) 

EeS^ti^Hi-^) V PN,E 

where Eq be the smallest energy shift such that E + Ea belongs to Spec(i/^*^^) for every E G Spec(i/(^)) and 7 is a 
suitable constant. The deterministic channel Cyes is given by the isometry 

V= Y. \M,E + Ea){N,E\. (D9) 

£;eSpec{ff(")) 

The fidelity is then given by 

F= Y Pm,e+e„> Y Pm.e+e„, ^min := min \E\. (DIO) 

BeSpec(i/(")) \E\<NE^ir, ^ 
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In the limit for large M, the shift Eq tends to zero and, by the central limit theorem, the r.h.s. converges to 
erf (iV_Eniin/ ^2A/ (if ^) ) , where erf denotes the error function. Lower-bounding the error function we then get 



F>l-\ ^ ^exp j^^\ (Dll) 



Let us now evaluate the probability of success, given by Pyes = 7 "^ EeSoedH'-'^)) P m ,e < 7 ■ Since tt^ in Eq. (D8) 
must be smaller than 1 for every E, the largest possible value for the constant 7 is 7 = niin^gspecfz/c™)) \ — ^^^^-^ — . 
In the hmit of large N, this gives Pyes < exp ( — iiH2)° ) 

3. Increasing the probability of success 

Here we devise super-replication processes that produce M = cN^^'^ copies for every desired e G [0, 1) and c > 0, 
with fidelity approaching 1 in the limit N — >■ 00, and with success probability lower bounded as Pyes > e^^ fW^ 
where f{N) is any fixed positive function such that limjv->.oo f{N) = 00. For e > 0, this proves in particular that 
one can achieve success probability larger than e~^ for every 6 > e. For e = 0, it proves in particular that one can 
achieve success probability larger than 1/N^ , for every (5 > 0. 

Consider the filter given by 



M,es- J2 7^e\N,E){N,E\, TTe 

BeSpec(H(")) 



Pn.e 

7 := mmiE\<E ' ^"'"^ 



^\ri^^ 1^1 ^ ^/w ■■= V^mwJjW) 



(D12) 



/(«) Y Pm, e+Eq 

\E\>Ef^r^). 



The combination of the filter Myes with the isometry V of Eq. (D9) gives a probabilistic process with fidelity 
lower bounded by F > J2\e\<e Pm,e+Eo^ which in the limit of large A^, the r.h.s. converges to the error function 

erf ( \/f{N)j , giving the bound 

i^>l-^47^exph/(iV)]. (D13) 

The probability of success of the above cloning process is lower bounded by 

Pyes> ^ TT%pN,E=l'^ ^ PM.E+Eg- (D14) 



In the hmit TV -> cx), the r.h.s. converges to gN,Ef(,^-,/gALEf^^-, erf f ^//(iV)j 



> g-^'/w. 



